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The Lagrangian density of Relativistic Hadrodynamics (RHD) is 
extended by introducing non- linear derivative (NLD) interactions 
of the nucleon with the meson fields. As a consequence, the nucleon 
selfenergy becomes both momentum and density dependent. With 
a single cut-off parameter, which regulates the NLD Lagrangian, 
our approach is compatible with results from microscopic nuclear 
matter calculations as well as with Dirac phenomenology. It also re- 
produces the correct behavior of the real part of the proton-nucleus 
potential both at low and at high proton energies. 
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1 Introduction 

Relativistic mean-field models have been widely accepted as a successful tool 
for the theoretical description of different nuclear systems such as nuclear mat- 
ter, finite nuclei and heavy-ion collisions [1]. The major advantage of relativis- 
tic mean-field (RMF) models has been a novel description of the saturation 
mechanism and an explanation of the strong spin-orbit force. An energy depen- 
dence of the Schrodinger equivalent optical potential [2] is naturally included 
as a consequence of a relativistic description. 

The common starting basis for RMF models is the relativistic hadrodynamics 
(RHD) with interacting hadrons, i.e., nucleons and mesons, as the relevant 
degrees of freedom. Starting from the pioneering work of Duerr [3] , simple La- 
grangians have been introduced [4], and since then many different extensions 
and treatments have evolved. The Lagrangian of RHD can be treated in dif- 
ferent approximations. For a practical application to nuclear matter and finite 
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nuclei one usually employs the mean-field approximation. In this approach the 
nucleon selfenergies become simple functions of various nucleon densities, but 
do not explicitly depend on energy or momentum of the nucleon. As a conse- 
quence a linear energy dependence of the optical potential arises. In principle, 
one can go beyond the mean-field approximation in a quantum field theoretical 
treatment by a systematic diagrammatic expansion. In this case, the nucleon 
selfenergies show a much more complicated structure that takes the non-trivial 
effects of the nuclear medium into account. In Dirac-Brueckner (DB) calcu- 
lations of nuclear matter based on a given nucleon-nucleon interaction it is 
possible to extract the corresponding nucleon selfenergies [5] . They depend on 
both the density of the nuclear medium and the energy and momentum of the 
nucleon. They reproduce the empirical saturation point of nuclear matter as 
well as the energy dependence of the optical potential for energies up to ~ 350 
MeV. 

However, for practical applications the DB approach is not suitable due to its 
high complexity and its apparent limitations at high energies and densities. 
Thus, as an alternative approach to these ab-initio strategies for the nuclear 
many-body system a phenomenological treatment of the problem in the spirit 
of the mean-field approximation has been established. The basic Lagrangian 
of RHD has to be modified for a quantitative description of static nuclear 
systems such as nuclear matter and/or finite nuclei. Therefore, it becomes 
mandatory to introduce new terms, e.g., including non-linear selfinteractions of 
the scalar [6] and vector [7] meson fields, or to modify existing contributions in 
the Lagrangian, e.g., introducing density dependent meson-nucleon couplings 
[8,9]. The model parameters have to be fitted to properties of nuclear matter 
and/or atomic nuclei, since they cannot be derived in a simple manner from 
a microscopic description. 

For a quantitative description of dynamical systems, e.g., hadron- and heavy- 
ion-induced reactions, the energy and momentum dependence of the nuclear 
mean-field becomes important. In conventional RMF models the meson-nucleon 
couplings are independent of the nucleon energy and momentum. This leads 
to a linear energy dependence of the Schrodinger equivalent optical potential 
with an unrealistically divergent behaviour at high energies. Indeed, analyses 
of proton-nucleus scattering data [2] show that the proton-nucleus optical po- 
tential levels off at proton energies of about 700 — 800 MeV. Thus, alternative 
RMF approaches have been developed by including additional non-local con- 
tributions, i.e., by introducing Fock-terms, on the level of the RMF selfenergies 
leading to a density and energy dependence [10]. However, such a treatment is 
not covariant and also its numerical realization in actual transport calculations 
is difficult [10]. Another approach has been proposed in [11] and more recently 
in [12] by introducing linear order derivative couplings in the Lagrangian of 
RHD. These additional terms lead to an energy (or momentum) dependence 
of a nucleon inside the nuclear matter, and modify also the nuclear matter 
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equation of state (EoS). Zimanyi and Moszkowski [11] studied such gradient 
terms with the conclusion of a softening of the nuclear EoS. A more detailed 
study by Typel [12] involved the investigation of both, the density dependence 
of the nuclear EoS and the energy dependence of the optical potential. While 
linear derivative interactions of meson fields with nucleons explain the em- 
pirical energy dependence of the optical potential, a soft EoS results from an 
introduction of an explicit density dependence of the nucleon-meson couplings 
with additional parameters. 

The purpose of the present work is to develop a covariant model, which gener- 
ates both, the correct density and momentum dependence of the selfenergies 
in an unified framework. Furthermore, we find that the density and the mo- 
mentum dependencies are indispensable and strongly correlated. The model 
proposed here is simple in realization and respects all the underlying symme- 
tries of the RHD Lagrangian. The approach followed here is to extend the 
original Lagrangian of RHD [4] by including non-Zmear derivative interactions 
of meson fields with nucleons. We refer it to as non-linear derivative (NLD) 
model, and it is presented in the Section 2. The NLD Lagrangian (Section 2) 
contains not only the nucleon spinor and its first derivative, but it depends on 
all higher order derivatives of the nucleon field. As a consequence, the stan- 
dard canonical formalism has to be generalized (Section 2.2). Field equations 
for nucleons and meson fields can be then derived (Section 2.3). The NLD 
model gives field equations with very simple structure in nuclear matter (Sec- 
tion 3) with density and energy dependent nucleon selfenergies. As a major 
advantage, both the equation of state (EoS) (density dependence) and the op- 
tical potential (energy dependence) are quantitatively well reproduced with a 
single parameter and the results are comparable with microscopic DB models, 
as in detailed discussed in Section 4. Finally, in Section 5, conclusions and an 
outlook complete the paper. 



2 General formalism 

This section describes the details of the formalism. In the first subsection the 
general Lagrangian density is constructed and discussed. In the second chapter 
various limiting cases arc investigated, before the Euler-Lagrange equations for 
the relevant degrees of freedom are formulated. 

2.1 The NLD Lagrangian 

The starting point is the Lagrangian density of RHD (RHD) [4], which de- 
scribes the interaction of nucleons through exchange of auxiliary meson fields 



3 



(Lorcntz-scalar, o", and the Lorentz- vector isoscalar and isovector meson fields 
uj^ and p^) in tiie spirit of tlie One-Boson-Excliange (OBE) model, see for in- 
stance [5] . In conventional RHD the interaction fields couple to the spinors via 
the corresponding Lorentz-density operators, e.g., ^^cr and (^^7^7^^^^^) 
for the scalar and vector sectors, respectively. In the mean-field approach to 
RHD the resulting mean-field potentials are only density dependent. For in- 
stance, its Lorentz- vector part rises linearly with nucleon density. 

In RHD the relevant degrees of freedom are nucleons and mesons (cr, a;, p) 
characterized by the corresponding free Lagrangians and an interaction part 
between them: 



+ T^^pP ^.P - -;G ^,yG 



The first line in Eq. (1) gives the (symmetrized) free Lagrangian for the nucleon 
field \[' = (^p, ^n)"^ with bare mass m. The other terms contain the standard 
Lagrangians for the scalar and vector mesons, cr, u)^ and p^, with the strength 
tensors F'^'^ = d'^u'' - d^u" and G^'' ^ d^p"" - d'^p^ for the isoscalar and 
isovector meson fields, respectively. 

We propose an extension of the standard RHD by introducing non-linear 
derivative couplings into the interaction Lagrangian density. Non-linear deriva- 
tives on the level of the Lagrangian are introduced by imposing an explicit de- 
pendence of the interaction Lagrangian density on higher order partial deriva- 
tives of the nucleon field ^: 



■'int 



(2) 



with the specified (symmetrized) interaction given by: 



^i: 



int ' 
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(3) 



The interaction between the spinor fields ^ , \E' and the meson fields has a 
similar functional form as in standard RHD [4]. However, now new operators 
T> acting on the nucleon fields appear, which is a generic non-linear function of 
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partial derivatives. The particular form of the non-linear derivative operator 
is not unique. Here a simple exponential form is assumed: 



The non-linear operators (4) explicitely depend on the partial derivate which 
is contracted with a dimensionless unit 4- vector v^. The choice of the vector 
field is not fixed from first principles. It can be identified with the nucleon 
4- velocity u^: 

= u^" := -^L= , (5) 



or with some other vector meson field: 



:= ^= . (6) 



No further specification is in this context necessary. We will indeed see that 
rearrangement contributions generated by in the corresponding field equa- 
tions cancel in nuclear matter. We thus consider it as an auxiliary unit vector. 
Furthermore, the additional factor e^^^ has been introduced for convenience 
so that one does not have to renormalize the conventional meson-nucleon cou- 
plings ga,uj,p of the RHD Lagrangian. The non-linear derivative operator T> 
contains A which will be interpreted as a cut-off parameter. Its value is sup- 
posed to be of natural hadronic scale of around 1 GeV. 

The entire interaction part of the Lagrangian density is symmetrized in order 
to keep the whole Lagrangian hermitian. It is invariant under global phase 
transformations which guarantees baryon number conservation. Translational 
invariance leads to a conserved energy-momentum tensor. 



2.2 Generalized formalism 



Before going forward with the field equations it is useful to study particular 
limiting cases of the general Lagrangian (3) for the nucleons. The isovector 
sector is not considered here for simplicity. A series expansion is here useful 
for studying limiting cases (i = cr, a;, p): 



5 



9ie A =9i2^ — 



n=o ^■ 

^g,- ^V'^l da datdf3±--- . 

The interaction terms of the Lagrangian density (3) in powers of the partial 
derivatives takes the form 

r _ /-(O) I I ^(2) I . . . ('7^ 

'"tnt '-int ~ '-int ~ '-int ~ ■ \' ) 



In zeroth order the interaction Lagrangian is 



Up to zeroth order in the expansion of the non-hnear derivative terms the 
standard Lagrangian of RHD is therefore retained: 
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-RHD 



(9) 



with minimal meson-nucleon coupling terms and introducing the covariant 
derivatives iD'^ := id^ — gi^uj^. Our general formalism contains as a limiting 
case, i.e., A — > oo, the conventional RHD. 

The consideration of higher orders leads to modified Lagrangians. E.g., up to 
first order a linear dependence on these derivatives is obtained, and so forth. 
However, a truncation of the expansion to an arbitrary order is not obvious. 
Instead, the most general case of a resummation of all higher non-linear orders 
in the interaction Lagrangian is meaningful, as introduced in Eq. (3). For this 
reason the full formalism will be considered for the rest of this work. 

The NLD Lagrangian >C is a functional of not only ^, ^ and their first deriva- 
tives, but it depends on all higher order space-time derivatives of the spinor 
fields ^ and ^. Thus the standard expressions for the Euler-Lagrange equa- 
tions and the Nocther Theorem do not apply. Instead of, we have to explic- 
itly derive the generalized Noether Theorem for a Lagrangian density which 
contains higher order derivatives, e.g., C = C{(j),da4>,dadp(t),- ■ ■). This step 
is important in defining conserved quantities such as the energy-momentum 
tensor and Noether currents. 

For a general functional containing fields and their higher order derivatives 
the Euler-Lagrange equations take the form: 
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For a general transformation of the fields and the coordinates the variational 
method for an arbitrary Lagrangian density (2), which contains higher order 
derivatives, leads to the following continuity equation (see appendix A for 
detailed derivations): 

^^[J^^-T^'-5x,]^Q. (11) 



where J'^ and Tf^^ are defined as: 



dC 



«9/3 + dpd^ 



+ 

+ 
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[9(9^0) 'd{d,dpct>) 
dC ^ dC 



dC ^ dC 
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d{di,dpd^(j)) 
±... 
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± 



dC ^ dC 



+ 
+ 



[9(9^0) 'd{d,dp<t>) 
dC ^ dC 



di,d^S(f)T 
dC 



d{d^df}d^(t)) 



_d{d^dj3(t)) d{d^d^d^(f)] 
dC ^ dC 



a 



diyd^d^d^ct)) d{d^d^df3d^(j)) 



±... 



(12) 



df3d^d''(f) + gi'''C . (13) 



In case of Poincare' transformations, i.e., S(j)T — (t>'{x') — (f){x) — and x' — 
x'^ + Sx^ we obtain d^T'^" = and therefore the conserved energy-momentum 
tensor density. In case of global phase transformations, i.e., Scpx ~ and 
Sx'^ = 0, the continuity equation d^J^^ = is retained. 

Contrary to the standard Noether theorem, now more terms proportional to 

higher order derivatives of the variation of the field S(j) appear, and for each 
order with respect to the derivative on 6(f), an infinite series appears. However, 
we will see that the infinite series can be resummed to compact expressions 
and that all the equations will simplify considerably. 
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2.3 Field equations for all degrees of freedom 



The application of the generahzed Euler-Lagrange equations (11) to the full 
Lagrangian density (3) with respect to the Spinor field ^ leads to a Dirac equa- 
tion with selfenergies, which in general contain infinite series with respect to 
the partial derivatives ^ . As shown in Appendix B, all series can be resummed 
to simple exponential functions leading to the following Dirac equation: 

[ji^iid^ - E'^) - (m - E,)] * = , (14) 



with Lorentz-vector and Lorentz-scalar selfenergies defined as: 

— v^i 8 Q +m —v(^ i 3 0-\-m 

E^ = g^uj'^e ^ + E(^ , E^ = g^ae x . (15) 



Both vector and scalar selfenergies E^, Eg show an explicit linear behavior 
with respect to the corresponding fields, cj^ and a, respectively, which looks 
like as in the linear Walecka model of RHD. It has been widely discussed in the 
hterature that the linear Walecka model does not reproduce the compressibil- 
ity of nuclear matter at saturation density, except if selfinteractions between 
the a meson fields [4,6] and the u fields [7] are introduced. One may thus 
expects the appearance of the same problem also in our formalism. On the 
other hand, the meson fields will contain a residual non-linear density behav- 
ior, which will be discussed below, together with the rearrangement term E(f. 
Apart from the density dependence, the selfenergies contain now a non-linear 
dependence on the partial derivatives, which in nuclear matter are related to 
the 4- momentum of a moving nucleon inside the nuclear matter at rest. 

Proca and Klein-Gordon equations for the meson fields are immediately de- 
rived from the Lagrangian density (3) by considering the explicit a;-dependence 
of the 4- vector (see again Eq. (6)): 



d.F'^'' + = -g^ 



+ ^^.^ (16) 



dad"'a + mla = -go 



-v'-'i a a+m 



(17) 



In order to keep the presentation of the model and the discussions transparent, 
we omit the p-meson contributions in all the expressions. It is, however, included in 
the calculations for neutron matter. 
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The Proca equation for the u meson field (16) shows a rather complex struc- 
ture of the source term indicating a highly non-linear behavior of the uo field 
(the rearrangement term is given in Eq. (B.8), however, it is not relevant 
as discussed below) . The same feature also holds for the source term of the a- 
meson field, although in the latter case its structure is simpler. Note that the 
meson field equations for the two mesons are coupled through the non-linear 
exponential terms. In nuclear matter these equations will simplify consider- 
ably, however, the coupling between both meson equations will remain. This 
feature together with the non-linearities in the source terms will generate a 
residual non-linear density behavior also for the vector field, which in the 
standard Walecka model is hnear in density. 

Depending on the particular choice of the the auxiliary unit vector in Eqs. (5) 
or (6) rearrangement terms appear in the vector part of the nucleon selfenergy 
(E(f) or in the source term of the Proca equation (il^). An important feature is, 
however, that independent of the choice of the rearrangement terms exactly 
cancel in nuclear matter (see appendix B.l and B.2) simplifying the equations 
of motion to a large extent. We, therefore, do not specify these rearrangement 
contributions here any further. 



3 Application to infinite nuclear matter 

Infinite nuclear matter is characterized by a homogeneous, isotropic and sta- 
tionary system at zero temperature. Due to these symmetries the field equa- 
tions simplify: the densities are constant in space-time, and only the 0-component 
of the currents survives. Consequently, the meson fields are constant numbers 
and the spatial components of the Lorentz-vector cj-field vanish. In the spirit 
of the mean-field approximation all meson fields are treated as classical fields 
implying the replacement of the source term operators in the meson field 
equations by their expectation values. Therefore = lo^ j ^JuJ^Jj^ = (1) 0) or 
^/i _ j ^j^ja = (1, 0), since the space-like components of the usual nucleon 
current j'^ —< ^7/^^ > vanish in average for nuclear matter at rest. 

In nuclear matter solutions of the spinor field are found with the plane wave 
ansatz (for simplicity isospin is neglected): 

*(s,p) = M(s,p)e-*^'''^'' , (18) 

with s and — (E.p) being the spin and 4-momentum of the nucleon, x'^ its 
time-space coordinate, and u the spinor field for positive energy states. 

In nuclear matter the rearrangement contributions cancel and the Dirac equa- 
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tion (14) maintains its form with selfenergies given by: 

E—m E—m 

= Quj^oe , S3 = g^ae a . (19) 

As one can see when the nuclcons arc at rest, E = m, the exponential factor 
is equal to unity and the equations are reduced to the ones from the Walecka 
model. The selfenergies show now a transparent behavior in energy. They con- 
tain a non-linear energy dependence, in particular, they exponentially decrease 
with increasing energy, as also expected from Dirac Phenomenology [2]. They 
are linear in the meson fields as in the linear Walecka model, however, this 
does not necessarily imply a linear density behavior of the mean field. 

Solutions of the Dirac equation are obtained by inserting (18) and (19) into 
the Dirac equation (14): 

7o£;** = (7-p + m*)*, (20) 



with 

E—m E—m . , 

E :— E — QujOJQe ^ , m :— m — QaO'e a . (^21) 



The Dirac equation (20) has the same structure as in the free space, however, 
with modified effective energy E* and effective (Dirac) mass m*. The spa- 
tial component of the 4-momentum p remains unchanged in nuclear matter. 
Solutions of the Dirac equation (20) are found in the usual way: 



/ 



u{s,p) = N 



(t>s 



(22) 



with spin eigenfunctions (ps and with the normalization N = y^-^^?-, which 
guarantees the orthogonality relation of the spinors. The following dispersion 
relation holds 

E*^-f = m*2 . (23) 



Note that, E* and m* depend both on density and, particular, on energy 
implying self consistency to obtain the dispersion relation E = E{p). 

In nuclear matter the NLD equations of motion for a and u simplify signifi- 
cantly to standard meson field equations: 

mluj^ = Qo^pQ , mla = gaps , (24) 
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with the following source terms {k — 4,2 for nuclear matter and neutron 
matter) : 



p,= (^e A ^) = __ j d?p—e A e(p-p^), 

Po = (*/e-^*) = ^ / rfV-^e(p - p^) . (25) 

Note that, the density po is not related to the conserved nucleon density pb- 
Later on, this quantity will be specified explicitely. Both vector and scalar 
densities po and Ps difi^er considerably from those of the conventional Walecka 
model due to the dispersion relation (23) and the appearance of the expo- 
nential terms, which act as damping factors. These effects lead to a coupled 
set of equations between the meson fields and the dispersion relation. They 
generate a residual non-linear density dependence for both the scalar field and 
especially the vector tu-meson. 

According to the generalized Noether-theorem (13) a conserved nucleon den- 
sity and an equation of state can be derived (see appendix C). The equation 

jO = P5 = (^/^) + ^(^70e-^^)a;o - ^(^e-^^)a (26) 



characterizes the conserved nucleon density p^, from which also the relation 
between the Fermi momentum pp and ps is uniquely determined. The expec- 
tation values po = (^7°e a~v1>) and ps = (^E'e ~'^) have been already 
evaluated in (25). The expectation value (\E'7°\E') = j^^i^n)^ is just the 
usual density of the Walecka model [4]. Unlike the Walecka model, which is 
recovered in the limiting case of A — > oo, the baryon density in the NLD ap- 
proach shows a more complex structure. This results from the application of 
the generalized Noether-theorem to the NLD Lagrangian in nuclear matter, 
in which the non-linear derivative terms are coupled with the meson fields, 
as shown in appendix C. The extra terms in Eq. (26) refiect the dressing of 
the nucleons with the meson clouds. The latter start to contribute at densities 
above saturation, e.g., at p^ = 2psat the nucleons contribute about 92% and 
the mesons 8% to the baryon density. 

The Equation of State (EoS) is obtained from the 00-component of the energy- 
momentum tensor T'^": 



T 



I -'■ / 2 2 2 2\ 



(27) 
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baryon density Pp [fm ] 



Fig. 1. Vector (oj, upper panel) and scalar (a, lower panel) meson fields as function 
of the nucleon density pB- Dashed: Linear Walecka model, solid: NLD model. 

with the additional expectation values: 



Ps 



K 



{2ny 



m 



{2ny 



Q Hi—m . , 



(28) 



4 Results 



The set of coupled equations (23,24) has to be solved self consistently together 
with the numerical integrations for the different densities (25,28). This proce- 
dure is done as function of the Fermi momentum pp and of the momentum p. 
The relation between the nucleon density pb and the Fermi momentum pf is 
given by Eq. (26). 

In principle, the NLD model contains no parameters. The original aNN and 
ujNN couplings can be taken from any linear Walecka model [4], as it has 
been done here. The cut-off parameter A has to be of natural size, i.e., of 
typical hadronic mass scale in this problem. In the following, A = 0.770 GeV 
is chosen for a quantitative discussion of the results. Of course, different cut-off 
parameters h.a,oj can be also used, which may be necessary when applying the 
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1000 
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Pb [fm"^] 

Fig. 2. Vector (upper panel) and scalar (lower panel) selfenergies as function of 
nucleon density. Dashed: Linear Walecka model, solid: NLD model, filled squares: 
DB model [5], filled circles: Sugahara and Toki [7]. 

NLD model to more complex nuclear systems such as finite nuclei. However, 
this is not the scope of the present work. 

We start the discussion with the dependence of the isoscalar meson fields 
as function of the nucleon density p^, see Fig. 1. In the standard linear 
Walecka model the a;-field shows the typical linear density behavior, while the 
scalar meson field saturates with increasing pB due to the suppression 7-factor 
m* / E* . In the NLD model, however, a saturation of the tu-field is observed 
too. This result is similar to that of Ref. [7], in which selfinteractions of the uj 
meson field were considered. This non-linear effect on the vector meson arises 
from the damping factor e~^/^, as discussed above. Non-linear effects to the 
scalar meson, on the other hand, are moderate due to partial compensation 
effects between the damping factor and the m* / E* term. 

The selfenerg proportional to the meson fields (see Eqs. (19)), how- 

ever, they are dominated by the factor e~^l^ which appears explicitly and is 
itself also density dependent through the in-medium dispersion relation. As 
can be seen in Fig. 2, the scalar part (lower panel) is thus further suppressed in 
the NLD model compared with the standard linear Walecka model. The den- 
sity dependence of the vector field in the NLD approach differs considerably 
from the usual linear dependence, in particular, at nucleon densities above 
saturation is strongly suppressed and tends to saturate at very high densi- 
ties. The agreement between the NLD model and microscopic DB calculations 
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Fig. 3. Equation of state for nuclear (left panel) and pure neutron matter (right 
panel). Dashed: Walecka model, solid: NLD model, filled squares: DB model [5]. 

[5] is remarkable over a wide region in nucleon density. Another possibility of 
suppressing the linear density dependence of the vector field was considered in 
Ref. [7], in which a selfinteraction term of the uj meson was introduced in the 
RMF Lagrangian. As can be seen again in Fig. 2, both models are very close 
to each other around saturation density and below, but differ considerably for 
PB > 0.3 fm-\ 

The non-linear density dependence significantly affects the equation of state 
(EoS), i.e., the binding energy per nucleon as function of nucleon density. This 
is demonstrated in Fig. 3 for nuclear and pure neutron matter. The standard 
linear Walecka model leads to an EoS with high stiffness resulting in a very 
high value for the compression modulus. The NLD model weakens the stiffness 
of the EoS for nuclear and pure neutron matter to a large extent. The agree- 
ment of the NLD-EoS with the underlying DB theory is successful indicating 
that saturation nuclear matter properties, e.g., binding energy per nucleon 
and compression modulus at saturation density, are fairly well reproduced by 
the NLD approach. It thus turns out that standard linear Walecka models 
without any explicit introduction of selfinteractions of the meson fields are 
able to describe nuclear matter properties by the introduction of non-linear 
derivatives with a single cut-off parameter. 

Similar effects are observed in the effective mass, as shown in the upper panel 
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Fig. 4. Upper panel: Nucleon density dependence of the effective (Dirac) mass m* 
in units of the bare nucleon mass M. Dashed: Walecka model, solid: NLD model, 
filled squares: DB model [5]. Lower panel: Energy dependence of the same quantity 
in the NLD approach. 

of Fig. 4. At saturation density a value of m* ~ 0.65m is obtained, which is 
again very close to DB predictions. The increase of the effective (Dirac) mass 
at high nucleon densities {ps > 0.3 fm~^) in the NLD model reflects the 
attractive character of the mean-field, which is consistent with the suppression 
of the scalar selfenergy at high densities and is quantitatively compatible with 
microscopic DB calculations. 

An important feature of the NLD model is the prediction of an energy (or mo- 
mentum) dependence of the nuclear mean-field with only one parameter. Dirac 
Phenomenology on elastic proton-nucleus scattering predicts a non-linear en- 
ergy dependence of the Schrodinger equivalent optical potential (see below), 
which cannot be reproduced in linear Walecka models nor in their extensions 
to non-linear meson field terms. The question arises if the NLD model can 
reproduce this feature with the same parameter A as used for the density 
dependence. 

In the following we consider the situation of a nucleon with particular momen- 
tum p (or kinetic energy -E^m) relative to nuclear matter at rest. The kinetic 
energy for an incident free nucleon with mass m and momentum p is usually 
defined as E^in = \/jP~+^rnF — m. In the nuclear medium Ekm reads [12,13] 




(29) 
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Fig. 5. NLD vector and scalar selfenergies (solid curves, as indicated) as function of 
the kinetic energy at fixed saturation density. The filled squares show DB calcula- 
tions taken from [13]. 

As a first impression Fig. 4 (panel on the bottom) shows the energy dependence 
of the effective (Dirac) mass at fixed nucleon density ps = Psat = 0.16 fm~^. 
In the standard linear Walecka model m* is independent on energy, since the 
scalar selfenergy is momentum independent. In the NLD model, however, the 
effective mass non-linearly rises with energy, and in the limit E oo (not 
shown in the figure) m* — >■ m is achieved, since e~^^^ 0. This result is in 
remarkable agreement with Dirac phenomenology [2]. 

The explicit energy dependence of the scalar and vector selfenergies in the 
NLD model are shown in Fig. 5 (solid curves). According to Eqs. (19) an 
approximately exponential decrease of Tis,v with increasing energy can be seen, 
which is unambiguously fixed by the same cut-off parameter A as used for 
the density dependence. The agreement with DB results is again remarkable 
in an energy region up to -Efc^ ~ 0.4 GeV. The DB model describes very 
well saturation properties of nuclear matter. It further gives a realistic energy 
dependence of the selfenergies in the considered energy region by using the 
parameters of the underlying bare iViV-potential [13]. The deviation of the 
fields between the NLD and DB approaches for energies above ~ 0.4 GeV will 
be important in describing the energy dependence of the optical potential. 

The energy dependence of the nuclear mean-field is usually determined by 
Dirac phenomenology in elastic nucleon-nucleus scattering [2]. The essential 
physical object here is the Schrodinger equivalent optical potential V^pt, which 
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Fig. 6. (left) Energy dependence of the Schrodinger equivalent optical potential 
at density psat = 0.16 fm~^. The theoretical calculations (dashed: linear Walecka 
model, solid: NLD model) are compared with data from Dirac phenomenology (filled 
squares) [2] and with DB calculations (filled stars) [13]. 

serves as a convenient means to characterize the in-medium interaction of a 
nucleon with momentum p relative to nuclear matter at rest. It is obtained 
by rearranging the Dirac equation (20) into a Schrodinger-like equation for 
the large (upper) component of the nucleon spinor, i.e., by performing a non- 
relativistic reduction of the Dirac equation. This procedure applied to Eq. (20) 
leads to a central Schrodinger-hke potential Vopt 



opt 



m 2m ^ ' 



(30) 



It is determined by the scalar and vector selfenergies and it rises linearly with 
energy, if the selfenergies do not depend explicitely on momentum. This is the 
case of the linear Walecka model, as can be seen in Fig. 6 (dashed curve). 
The DB model (filled stars), on the other hand, reproduces the empirical 
behavior of the optical potential at low energies, since the parameters of the 
underlying free NN- interaction are fitted to low energy scattering data [13]. 
The NLD model (solid curve) with its non-linear energy dependence weakens 
the linear stiffness of the original Walecka model to a large extent, and the 
empirical energy behaviour of the optical potential can be reproduced without 
the introduction of any further parameters. 

The strong suppression of the mean field at higher energies at saturation 
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Fig. 7. Energy dependence of the Schrodinger equivalent optical potential at dif- 
ferent densities, as indicated. For each density both calculations are shown, linear 
Walecka model (dashed) and NLD approach (solid). The potentials at ps = O.bpsat 
(dotted lines), pB = Psat (solid lines) and pB = 2psat (dot-dashed lines) are shifted 
downwards by 100 MeV, 50 MeV and MeV, respectively. Horizontal lines mark 
the range in the energy where the potentials cross zero. 

density is maintained also at higher nucleon densities, as shown in Fig. 7. For 
highly compressed matter, as it occurs in high energy heavy ion collisions, one 
thus expects a significant suppression of the repulsion of the nuclear mean field 
in the NLD approach. The novel result of the NLD model is thus a softening of 
the EoS at high densities and the Schrodinger equivalent optical potential at 
high energies and also high compressions. These features will be of particular 
interest for the theoretical description of high energy heavy-ion collisions, as 
they will be investigated at the compressed baryonic matter experiment at the 
new FAIR facihty at GSI [14]. 



5 Conclusions and outlook 



In this work the conventional linear Walecka model of relativistic hadrody- 
namics was extended by introducing non-linear derivative operators in the 
interaction Lagrangian. A generalization of standard field theoretical tech- 
nics for Lagrangians depending on all higher order derivatives of a field was 
necessary in deriving field equations and applying the Noether theorem in a 
consistent way. Starting from rather complicated field equations it was possi- 
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blc to reduce them to compact and simple expressions by resuming all higher 
order derivatives. 



The model was applied in nuclear matter using a single parameter of natural 
hadronic mass scale. Astonishing were the non-hnear density effects on the 
vector field without the introduction of any additional selfinteraction terms 
in the original Lagrangian of the linear Walecka model. These effects lead to 
a softening of the equation of state for nuclear and pure neutron matter at 
high densities. It was possible to describe quantitatively well the empirically 
known saturation properties. The results were also comparable to predictions 
of microscopic DB calculations over a wide density range. The approach lead 
furthermore to a momentum dependence of the selfenergies. Most noticeable 
was the correct reproduction of the energy dependence of the Schrodinger 
equivalent optical potential by utilizing the same parameter. The microscopic 
origin or the underlying mechanism of the desired results is not relevant at 
the present level of study. 

In principle, one could go beyond the present study and apply the model to 
more complex nuclear systems, e.g., finite nuclei. In this context the present 
model can be further extended by introducing different cut-off parameters 
or different non-linear operator functions than the exponential form used in 
this work. Another physical object of interest would be heavy ion collisions 
at high energies with highest baryon compression, as they are planned at the 
new FAIR facility at GSI, in which the density and, in particular, the energy 
dependence of the nuclear mean-field play an important role. The applica- 
tion of the present model to heavy-ion collisions in the spirit of a covariant 
transport theory based on the present Lagrangian would be a great challenge 
for the future in studying hadronic matter under extreme conditions with the 
ultimative goal of exploring the equation of state at supra-normal densities. 

Acknowledgments: This work is supported by BMBF. 

Appendix 

A Generalized Euler-Lagrange equations and Noether theorem 

The procedure to derive the generalized Euler-Lagrange equations follows the 
canonical formalism and starts from the variation of the action. Standard steps 
like the integration by parts lead to the following expression: 




a 



dC 



d{dadi3(p) 



dC 



H h 



(A.l) 



= . 
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The derivation of the conserved energy-momentum tensor and current is more 
involved. The conserved energy-momentum tensor and current are obtained 
from a variation of the Lagrangian with respect to infinitesimal changes of the 
field (j) and its argument. The total variation of the field reads 

St(P{x) := (P'{x') - (P{x) - - + - 0(x)] (A.2) 



which in first order reduces to the usual form 

5T(t>{x) — 54>{x) + da(p{x)5x°' 



(A.3) 



with 5(f) — (j)' — 4> and 5x^ being the small variation of the field at fixed 
argument and the infinitesimal change of it, respectively. 

The variation of the Lagrangian density functional thus reads 



C[ct)'{x'), d^dpct^'ix'), ■■■,]- C[ct){x'), d^dpct){x'), • • • , ] 

+C[<p{x'), da(j){x'), dadf3(j){x'), •••,]- C[(j){x), da(j){x), dadj3(t){x), •••,] = . 

(A.4) 

The first line in Eq. (A.4) is just the variation of C with respect to the fields 
0(3;), da<p- ■ ■ • at fixed argument, whereas the second line gives the variation 
of the Lagrangian with respect to the argument which can be evaluated up to 
first order. In total we obtain: 



5C + d^CSx"" - 



(A.5) 



with 



5C 



dC 



dC dC 



d{d^(t>) 



d{dadf3(f)) 



+ 



(A.6) 



Inserting in the first term in Eq. (A.6) the Euler- Lagrange equations (11), we 
obtain after the application of the product rule following expression: 



■9/3 + dpd^ 



dC 



+ 



+ 



dC 
d{d^dy(j)) 
dC 



a 



d{df,dyd^(, 
+ • • •} . 



d{d^dp(p) 
dC 

"'d{d^d0d,^\ 
a 



d{d^dpd^(p) 



±... 



(A.7) 



20 



As next step we insert the expression A. 7 into (A. 5) obtaining: 

da (J" - T'"'6x^) = 



(A.8) 



which is the desired result, see also (11). 



B Field equations in the NLD model 



In this appendix the most important steps leading to the field equations for 
the spinor ^' and the meson fields a and uu are given. The starting basis for 
the derivation of a Dirac equation from the Lagrangian density (3) are the 
generahzed Euler-Lagrange equations (11). We give the terms up to second 
order explicitly. 

In order to show the disappearance of the rearrangement terms in nuclear 
matter, we choose in the derivation of the Dirac equation for the auxiliary 
unit vector the usual nucleon current, see Eq. (5), in which rearrangement 
terms only in the vector selfenergy will appear. For the derivation of the Proca 
equation {u meson) v'^ = lo^ j \J uj qUJ'^ will be chosen, in order to show that 
the dissapereance of the rearrangement terms in nuclear matter is of general 
nature, i.e., independent of the particular choice of the auxiliary vector . 



B.l Dirac equation 



For the derivation of the Dirac equation wc start from the generalized Euler- 
Lagrange equations (A.l). We derive explicitely the first three orders in the 
derivative expansion: 



dc dcdv^' 1 -OM, 

-= + -= = -^aid^ib - m'^ 



1 



2^. 



*(T + e A 



i 3 f-iV' +m 



( ^(^ 9Je^^7"^a;„ + *7"e- 



-i d^)'^U}a 



X 



7^* j^juY'^ 
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dC 



dC 



2 ' ^ 2 A ' 2 A 



2 (a)^ 2! 



2(A)2 2! 



(B.l) 



Inserting (B.l) into (A.l) we obtain after some trivial algebra (up to deriva- 
tives of the meson fields which in any case vanish in nuclear matter) : 



1 



A 



1 1 
2! A2 
1 

2!A2 



^ Vaid'^ + ^^v^vpid'^id^ ip . . + E^^ = . (B.2) 



with the rearrangement contribution: 



YX I 9^)e^^7'*^5 + *7'e ^ (-i 9^)*u;5 



X 



7^ 



^3 



(B.3) 



Note that, the derivatives inside the brackets act to the spinors on the right. 
The terms inside the brackets represent Taylor expansions which can be re- 



sumed to the exponential function e 
Dirac equation 



A . This leads us to the following 
[j^iid^^ - E'^) - (m - E,)] * = , (B.4) 



as it was written in Eq. (14). 

In infinite nuclear matter only the 0-component of the uu meson field survives, 
and in the mean-field approximation all density operators are replaced by 
expectation values. Furthermore only the 0-component of the usual nucleon 
current j'^ — (^'7'*^) remains. Thus, the rearrangement contribution reads: 



,(B.5) 
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and by evaluating the contraction between the expressions inside the two 
brackets we obtain: 



Qi^ n T \ 9o" E-m 



(B.6) 



The first term in Eq. (B.6) is trivially zero. In the second line of Eq. (B.6) 
both expectation values are odd functions of the vector p and therefore they 
vanish. 



B.2 Meson field equations of motion 



For the meson fields a and uj the derivations are straightforward, since here 
one has to use the standard Euler- Lagrange equations. The following Proca 
equation which includes the rearrangement and the source terms is readily 
derived: 



A + ^T^e A ^ 



+ (B.7) 



The rearrangement term in the Proca equation appears if for the auxiarily 
vector = cu^/ ^JlOoLo'^ is chosen. It takes the following form: 



:= - 



1^ 
2 A 



i 8 av'-'+m 



id 



id UaOj" 



id 



+ 



ijjnijj" id 



1^ 

2 A 



•q/e A- 



id 



id uOnUO^ 



id 



+ 



uJni^'^ id 



(B.8) 



In nuclear matter and in the mean-field approximation the rearrangement 
contribution fi^ to the vector meson vanishes for the same reasons as for that 
rearrangement term for the nucleon field (see previous subsection). 
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C Derivation of the Noether current 



Here the most important steps leading to the Noether current (26) starting 
from the general expression (12) are explicitly shown. The derivation of the 
encrgy-momcntTim tensor (13) is then straightforward. We show that all higher 
order derivative terms can be resumed to simple exponential expressions for 
nuclear matter at rest in the mean-field approximation. 

The derivative of the Lagrangian density (3) with respect to higher order par- 
tial derivatives has been demonstrated in the appendix B.l, see in particular 
Eq. (B.l), thus we don't repeat this intermediate step. 

We consider the Noether current in the NLD model starting from the general 
expression (12) and taking only terms up to third order. With — —ie^ 
and — ie^ we obtain: 



ie 



+ 
+ 



dC dC 

a(^ ^d{d,d^^) 

dC dC 

dC DC 



8C 



8a^ 



8C 



-8, 



8C 



8{8,^) "^'8{8,8p^>) ' 8{8,8p8,^) 
dC „ dC 



8C 



+ 8p^ 
+ 8,8^^> 



d{d,dp^) ^'8{8,8pd,^)\ 

dC ^ dC 

— o. 



(C.l) 



We rewrite now Eq. (C.l) by separating the terms between the different orders 
in the partial derivatives: 



le 



0(0) + 0(1) + 0(2) ^ . . . ^ 0(n) 



(C.2) 



with the different contributions given by: 



0(0) = ( _ ^. 



(C.3) 
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dC 



+ 



dc • 

d{d^da^) 
__dC__ 



(C.4) 



'd{d,d^d,^>) ' ' d{d,dpd,^) ) 

( dC, - dC \ 



In mean-field approximation it turns out that first order terms (in energy) 
appear twice, second order terms appear three times, and so forth: 

0(0) = (fr/^) + ^(fr/^)^Q _ ^{^>m)a, (C.6) 

O(^) = f|(^/|^^)a;o-||(^|^*)a, (C.8) 



and for the n — th-order one obtains 

^(n-l) ^ (-)"^^(^/f^*)-0 - (-)'^^^(^S*>- ■ (C.9) 



An!' A"-i ' ' ' A n! ' A'^ 

Therefore all the infinite series of terms can be resumed into simple expres- 
sions. Adding all these terms and considering the fact that ^ = -r^-rrr we 

° ° n! {n—l)' 

obtain the desired result, see Eq. (26). 
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